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Within this paper, control-relevant vehicle design concepts are examined using a widely 
used 3 DOF (plus flexibility) nonlinear model for the longitudinal dynamics of a generic 
carrot-shaped scramjet powered hypersonic vehicle. The impact of elevator size and place- 
ment on control-relevant static properties (e.g. level-flight trimmable region, trim controls, 

AOA, thrust margin) and dynamic properties (e.g. instability and right half plane zero as- 
sociated with flight path angle) are examined. Elevator usage has been examine for a class 
of typical hypersonic trajectories. 

I. INTRODUCTION AND OVERVIEW 

Motivation. With the historic 2004 scramjet-powered Mach 7 and 10 flights of the X-43A 1-4 , hypersonics 
research has seen a resurgence. This is attributable to the fact that air-breathing hypersonic propulsion is 
viewed as the next critical step toward achieving (1) reliable, affordable, routine access to space, as well 
as (2) global reach vehicles. Both of these objectives have commercial as well as military implications. 
While rocket-based (combined cycle) propulsion systems 5 are needed to reach orbital speeds, they are much 
more expensive to operate because they must carry oxygen. This is particularly costly when traveling at 
lower altitudes through the troposphere (i.e. below 36,152 ft). Current rocket-based systems also do not 
exhibit the desired levels of reliability and flexibility (e.g. airplane like takeoff and landing options). For 
this reason, much emphasis has been placed on two-stage-to-orbit (TSTO) designs that involve a turbo- 
ram-scramjet combined cycle first stage and a rocket second stage. This paper focuses on control challenges 
associated with scramjet-powered hypersonic vehicles. Such vehicles are characterized by significant aero- 
thermo-elastic-propulsion interactions and uncertainty 1-17 . 

Controls-Relevant Hypersonic Vehicle Modeling. The following significant body of work (2005- 
2009) 7-9, 18-27 examines aero-thermo-elastic-propulsion modeling and control issues using a first principles 
nonlinear 3-DOF longitudinal dynamical model which exploits inviscid compressible oblique shock-expansion 
theory to determine aerodynamic forces and moments, a ID Rayleigh flow scramjet propulsion model with 
a variable geometry inlet, and an Euler-Bernoulli beam based flexible model. The vehicle is 100 ft long with 
weight (density) 6232 lb per foot of depth and has a bending mode at about 20 rad/sec. The controls include: 
elevator, stoichiometrically normalized fuel equivalency ratio (FER), diffuser area ratio (not considered in 
this work), and a canard (not considered in this work). A more complete description of the vehicle model 
can be found in previous works 7,28 . 

More recent modeling efforts have focused on improved propulsion modeling 29,30 that captures precombus- 
tion shocks, dissociation, wall heat transfer, skin friction, fuel-air mixing submodel, and finite-rate chemistry. 
The computational time associated with the enhanced model is significant, thus making it cumbersome for 
control-relevant analysis. The simple ID Rayleigh flow engine model discussed within 7, 18,25,28 will be used 
in the current paper. 
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TSTO Vehicle Trajectory. TSTO launch vehicles consist of a turbo-ram-scramjet combined cycle first 
stage launch vehicle and a rocket powered orbiter. 31 A combination of turbo-rocket-ram powered engines 
propel the craft from the ground to roughly Mach 5. An air-breathing scram-jet engine accelerates the vehicle 
to Mach 8 along a dynamic pressure of 2000 psf. High dynamic pressures result in better performance for the 
air-breathing engine. Structural limitations require staging (separation of the launch and orbiter vehicles) 
to occur at a much lower dynamic pressure, roughly 200 psf. An pull-up maneuver is performed at Mach 
8 to increase the altitude until the desired dynamic pressure. After staging, the orbiter vehicle is propelled 
to low earth orbit by a rocket engine, while the launch vehicle begins a un-powered descent. This work will 
focus on the altitude pull-up maneuver. 

Hypersonic Vehicle Control Issues. Within this paper, we exploit the generic carrot-shaped vehicle 
3DOF (plus flexibility) model presented in [7,18,25,28,32,33]. A myriad of issues exist that make control 
design for this hypersonic vehicle a potentially challenging problem: 

• Input/Output Coupling. For this system, velocity control is achieved via the FER input. Flight path 
angle (FPA) control is achieved with the elevator 34 . However, there is significant coupling between 
FER and FPA. 

• Unstable/Nonminimum Phase. Tail controlled vehicles are characterized by a non-minimum phase 
(right half plane, RHP) zero that is associated with the elevator to FPA map 27 . This RHP zero limits 
the achievable elevator-FPA bandwidth (THE) 35-37 . In addition, the rearward situated scramjet and eg 
(center of gravity), implies an inherent pitch-up vehicle instability. This instability requires a minimum 
BW for stabilization 28 . To address these potentially conflicting specifications, one approach has been 
to exploit the addition of a canard 18,34,38-40 . It is understood, of course, that any canard approach 
would face severe heating, structural, and reliability issues. 

• Varying Dynamic Characteristics. Within 28 , it is shown that the nonlinear model changes significantly 
as a function of the flight condition. Specifically, it is shown that the vehicle pitch-up instability and 
non-minimum phase (NMP) zero vary significantly across the vehicle’s trimmable region. In addition, 
the mass of the vehicle can be varied during a simulation in order to represent fuel consumption. 
Several methods have been presented in the literature to deal with the nonlinear nature of the model. 

Papers addressing modeling issues include: nonlinear modeling of longitudinal dynamics 27 , heating 
effects and flexible dynamics 9, 23, 41 , FPA dynamics 38 , unsteady and viscous effects 8, 19 , and high fidelity 
engine modeling 29,30,42 . 

Papers addressing nonlinear control issues include: control via classic inner-outer loop architecture 43 , 
nonlinear robust /adaptive control 34 , robust linear output feedback 40 , control-oriented modeling 18 , lin- 
ear parameter- varying control of flexible dynamics 44 , saturation prevention 21, 45, 46 , and thermal choking 
prevention 28,46 . 

• Uncertain Flexible Modes and Coupling to Propulsion. Flexible dynamics have been captured within 
the model by approximating a free- free Euler-Bernoulli beam using the assumed modes method 23 . 
Three flexible modes are used to approximate the structural dynamics. A damping factor of £ = 0.02 
is assumed. The associated mode frequencies are ui i = 20.94 rad/sec, W2 = 50.58 rad/sec, W3 = 100.5 
rad/sec. These modes must be adequately addressed within the control system design process. While 
performance can be improved by increasing controller complexity (e.g. higher order notches) 44 , one 
must be wary of, and careful in dealing with, modal/damping uncertainty issues. This is particularly 
important because structural flexing impacts the bow shock. This, in turn impacts the scramjet’s inlet 
properties, thrust generated, aft body forces, the associated pitching moments, and hence the vehicle’s 
attitude. Given the tight altitude-Mach flight regime - within the air-breathing corridor 5 - that such 
vehicle must operate within, the concern is amplified. In short, one must be careful that the control 
system BW and complexity are properly balanced so that these lightly damped flexible modes are not 
overly excited. 

• Control Saturation Constraints. Control saturation is of particular concern for unstable vehicles such 
as the one under consideration. Two specific saturation nonlinearities are a concern for any control 
system implementation. 


2 of 23 


American Institute of Aeronautics and Astronautics 



— Maximum Elevator/Canard Deflection and Instability. FPA is controlled via the elevator/canard 
combination 38 . Because these dynamics are inherently unstable, elevator saturation can result in 
instability 45 . Classical anti- windup methods may be inadequate to address the open loop unstable 
vehicle. The constraint enforcement method within 45,47 and generalized predictive control 48 have 
been used to address such issues. This paper examines several methods to prevent instability due 
to elevator saturation. It should be noted that control surface/actuator rate limits must also be 
properly addressed by the control system in order to avoid instability. 

— Thermal Choking/Unity FER: State Dependent Constraint. As heat is added within the combus- 
tor, the supersonic air flow is slowed. If enough heat is added, the combustor exit Mach number 
will approach unity, and the flow is said to be thermally choked 49 . If additional heat is added, 
the upstream conditions can be altered. This can (in principle) lead to engine unstart 5 - a highly 
undesirable condition. The amount of FER that causes thermal choking at a particular flight 
condition is referred to as the thermal choking FER, or FERtc- In general, FERtc depends 
upon the free-stream Mach, free-stream temperature, pressure, and density (which depend on the 
altitude), and the flow turn angle (vehicle geometry + AOA + elastic deflection) 28,48 . In addition, 
since the model does not capture what happens when FER > l 27 , it is natural to restrict FER 
below unity. Given the above, it follows that the minimum of these two constraints dictates the 
available FER at a given flight condition. The resulting state dependent FER constraint can be 
computed (on-line) based on the flight condition, and must be accounted for by the control law. 
Here, uncertainty is of great concern because of the potential unstart issues - issues not captured 
within the model. Engineers, of course, would try to “build-in protection” so that this is avoided. 
As such, engineers are forced to tradeoff operational envelop for enhance unstart protection. 

Control-Relevant Vehicle Design Issues. Despite the successful integrated approach taken by the X- 
43A team, as well as other prior successful flight control efforts, far too often aerospace vehicle design has not 
significantly involved the discipline of controls until very late in the vehicle design process or even afterwards. 
Research programs over the past two decades have suggested that for the anticipated hypersonic vehicles, the 
traditional “sequential” approach is not likely to work. This is attributable, in part, to complex uncertain 
nonlinear coupled unstable, non-minimum phase, flexible dynamics together with stringent flight corridor 
and variable constraints (e.g. specific impulse, fuel use, maximum dynamic pressure, engine temperatures 
and pressures). For such vehicles, an integrated multidisciplinary “parallel” approach - involving multiple 
disciplines up front - is essential. This is particularly true when tight flight control specifications must be 
satisfied in the presence of significant uncertainty. 

Goals and Contributions of Paper. This paper addresses a elevator surface issues that are of concern 
to both vehicle and control system designers. In short, this paper represents a step toward answering the 
following critical control-relevant vehicle design questions: 

1. How do elevator design properties impact a vehicles static and dynamic properties? 

2. How do these impact control system design? 

3. How should the elevator be designed in the presence of known elevator saturation? 

In short, this paper illustrates fundamental tradeoffs that vehicle and control system designers should jointly 
consider during the early stages of vehicle conceptualization/design. While vehicle designers may want to use 
a higher fidelity model (e.g. Euler based CFD with boundary layer reconstruction or Navier-Stokes based 
CFD50) to conduct more accurate vehicle trade studies, this paper shows that a (first principles) 3DOF 
nonlinear engineering model - such as that used in the paper - may be very useful during the early stages of 
vehicle conceptualization and design. 

Organization of Paper. The remainder of the paper is organized as follows. 

• Section II provides an overview of the dynamical model to be used in our studies. 

• Section III demonstrates the tradeoffs for elevator surface sizing and placement. 

• Section IV outlines the control and guidance architectures used to perform the pull up maneuver. 

• Section IV. A contains nonlinear simulations of the pull up maneuver. 
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• Section V summarizes the paper and presents directions for future research. 


II. DESCRIPTION OF NONLINEAR MODEL 

In this paper, we consider a first principles nonlinear 3-DOF dynamical model for the longitudinal dy- 
namics of a generic scramjet-powered hypersonic vehicle '~ 9, 18-2 ' . The vehicle is 100 ft long with weight 
(density) 6,154 lb per foot of depth and has a bending mode at about 22 rad/sec. The controls include: 
elevator, stoichiometrically normalized fuel equivalency ratio (FER), diffuser area ratio (not considered in 
our work), and a canard (not considered in our work). The vehicle may be visualized as shown in Figure l 8 . 



Modeling Approach. The following summarizes the modeling approach that has been used. 

• Aerodynamics. Pressure distributions are computed using inviscid compressible oblique-shock and 
Prandtl- Meyer expansion theory 10, 15, 27, 49 . Air is assumed to be calorically perfect; i.e. constant specific 
heats and specific heat ratio 7 = f ^ = 1.4 10,49 . A standard atmosphere is used. 

Viscous drag effects (i.e. an analytical skin friction model) are captured using Eckerts temperature 
reference method 8,10 . This relies on using the incompressible turbulent skin friction coefficient formula 
for a flat plate. Of central importance to this method is the so-called wall temperature used. The 
model assumes a nominal wall temperature of 2500°R 8 . 

Unsteady effects (e.g. due to rotation and flexing) are captured using linear piston theory 8,50 . The idea 
here is that flow velocities induce pressures just as the pressure exerted by a piston on a fluid induces 
a velocity. 

• Propulsion. A single (long) forebody compression ramp provides conditions to the rear-shifted scramjet 
inlet. The inlet is a variable geometry inlet (variable geometry is not exploited in our work). 

The model assumes the presence of an (infinitely fast) cowl door which uses AOA to achieve shock- 
on-lip conditions (assuming no forebody flexing). Forebody flexing, however, results in air mass flow 
spillage 2 '. At the design cruise condition, the bow shock impinges on the engine inlet (assuming no 
flexing). At speeds below the design-flight condition and/or larger flow turning angles, the cowl moves 
forward to capture the shock. At larger speeds and/or smaller flow turning angles, the bow shock is 
swallowed by the engine. In either case, there is a shock reflected from the cowl or within the inlet 
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(i.e. we have a bow shock reflection). This reflected shock further slows down the flow and steers it 
into the engine. It should be noted that shock-shock interactions are not modeled. 

The model uses liquid hydrogen (LH2) as the fuel. It is assumed that fuel mass flow is negligible 
compared to the air mass flow. Thrust is linearly related to FER for all expected FER values. For 
large FER values, the thrust levels off. In practice, when FER > 1, the result is decreased thrust. This 
phenomena 27 is not captured in the model. As such, control designs based on this nonlinear model (or 
derived linear models) should try to maintain FER below unity. 

The model also captures thermal choking. A (state dependent) saturation level - associated with FER 
(e.g. thermal choking and unity FER) - and a useful FER margin definition (one that can be used for 
the design of control systems for scramjet-powered hypersonic vehicles) was addressed in 43 . Finally, it 
should be noted that the model offers the capability for addressing linear fuel depletion. 

• Structural. A single free-free Euler-Bernoulli beam partial differential equation (infinite dimensional 
pde) model is used to capture vehicle elasticity. As such, out-of-plane loading, torsion, and Timoshenko 
effects are neglected. The assumed modes method (based on a global basis) is used to obtain natural 
frequencies, mode shapes, and finite-dimensional approximants. This results in a model whereby the 
rigid body dynamics influence the flexible dynamics through generalized forces. Within the current 
model, forebody deflections influence the rigid body dynamics via the bow shock which influences 
engine inlet conditions, thrust, lift, drag, and moment 23 . Aftbody deflections influence the AOA seen 
by the elevator. As such, flexible modes influence the rigid body dynamics. 

The nominal vehicle is 100 ft long. The associated beam model is assumed to be made of titanium. 
It is 100 ft long, 9.6 inches high, and 1 ft wide (deep). This results in the nominal modal frequencies 
u>i = 21.02 rad/sec, u>2 = 50.87 rad/sec, L 03 = 101 rad/sec. When the height is reduced to 6 inches, 
then we obtain the following reduced modal frequencies: uq = 10.38 rad/sec, uq = 25.13 rad/sec, 
W3 = 49.89 rad/sec. Future work will examine vehicle mass-flexibility-control trade studies 18 . 

• Actuator Dynamics. Simple first order actuator models (contained within the original model) were 

used in each of the control channels: elevator - — , FER - s+ ° 0 , canard - (Note: canard not 

used in our study). An elevator saturation of ±30° was used. 21,45 Elevator position and rate saturation 
become very important given the vehicle’s (open loop) unstable dynamics. 

Generally speaking, the vehicle exhibits unstable non-minimum phase dynamics with nonlinear aero-elastic- 
propulsion coupling and critical (state dependent) FER constraints. The model contains 11 states: 5 rigid 
body states (speed, pitch, pitch rate, AOA, altitude) and 6 flexible states. 

Unmodeled Phenomena/Effects. All models possess fundamental limitations. Realizing model limita- 
tions is crucial in order to avoid model misuse. Given this, we now provide a (somewhat lengthy) list of 
phenomena/effects that are not captured within the above nonlinear model. (For reference purposes, flow 
physics effects and modeling requirements for the X-43A are summarized within [51].) 

• Dynamics. The above model does not capture longitudinal-lateral coupling and dynamics 52 and the 
associated 6DOF effects. 

• Aerodynamics. Aerodynamic phenomena/effects not captured in the model include the following: 
boundary layer growth, displacement thickness, viscous interaction, entropy and vorticity effects, lam- 
inar versus turbulent flow, flow separation, high temperature and real gas effects (e.g. caloric im- 
perfection, electronic excitation, thermal imperfection, chemical reactions such as O2 dissociation) 10 , 
non-standard atmosphere (e.g. troposphere, stratosphere), unsteady atmospheric effects 6 , 3D effects, 
aerodynamic load limits. 

• Propulsion. Propulsion phenomena/effects not captured in the model include the following: cowl 
door dynamics, multiple forebody compression ramps (e.g. three on X-43A 53,54 ), forebody boundary 
layer transition and turbulent flow to inlet 53,54 , diffuser losses, shock interactions, internal shock effects, 
diffuser-combustor interactions, fuel injection and mixing, flame holding, engine ignition via pyrophoric 
silane 3 (requires finite-rate chemistry; cannot be predicted via equilibrium methods 55 , finite-rate chem- 
istry and the associated thrust- AO A-Mach-FER sensitivity effects 30 , internal and external nozzle losses, 
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thermal choking induced phenomena (2D and 3D) and unstart, exhaust plume characteristics, cowl 
door dynamics, combined cycle issues 5 . 

Within [30], a higher fidelity propulsion model is presented which addresses internal shock effects, 
diffuser-combustor interaction, finite-rate chemistry and the associated thrust- AO A-Mach-FER sensi- 
tivity effects. While the nominal Rayleigh-based model (considered here) exhibits increasing thrust- 
AOA sensitivity with increasing AOA, the more complex model in 30 exhibits reduced thrust-AOA 
sensitivity with increasing AOA - a behavior attributed to finite-chemistry effects. 

Future work will examine the impact of internal engine losses, high temperature gas effects, and 
nozzle/plume issues. 

• Structures. Structural phenomena/effects not captured in the model include the following: out of plane 
and torsional effects, internal structural layout, unsteady thermo-elastic heating effects, aerodynamic 
heating due to shock impingement, distinct material properties, 56 and aero-servo-elasticity 57,58 . 

— Heating-Flexibility Issues. Finally, it should be noted that Bolender and Doman have addressed a 
variety of effects in their publications. For example, within [9,23] the authors address the impact of 
heating on (longitudinal) structural mode frequencies and mode shapes. Comprehensive heating- 
mass-flexibility-control studies will be examined further in a subsequent publication. 

• Actuator Dynamics. Future work will examine the impact of actuators that are rate limited; e.g. ele- 
vator, fuel pump. 

It should be emphasized that the above list is only a partial list. If one needs fidelity at high Mach numbers, 
then many other phenomena become important; e.g. O 2 dissociation 10 . 


Longitudinal Dynamics. The equations of motion for the 3DOF flexible vehicle are given as follows: 
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where L denotes lift, T denotes engine thrust, D denotes drag, A4 is the pitching moment, N. i: denotes 
generalized forces, Q demotes flexible mode damping factor, 07 denotes flexible mode undamped natural 
frequencies, m denotes the vehicle’s total mass, I yy is the pitch axis moment of inertia, go is the acceleration 
due to gravity at sea level, and Re is the radius of the Earth. 


• States. Vehicle states include: velocity v, FPA 7, altitude h, pitch rate q, pitch angle 9 , and the flexible 
body states 7?i , rji , 72 , 72 , 7 , 3 , 7.3 ■ These eleven (11) states are summarized in Table 1. 

• Controls. The vehicle has three (3) control inputs: a rearward situated elevator S ei a forward situated 
canard S c a , and stoichiometrically normalized fuel equivalence ratio (FER). These control inputs are 
summarized in Table 2. In this paper, we will only consider elevator and FER; i.e. the canard has been 
removed. 


a In this paper, we have removed the canard. Future work will examine the potential utility of a canard as well as its viability. 
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« 

Symbol 

Description 

Units 

1 

V 

speed 

kft/sec 

2 

7 

flight path angle 

deg 

3 

a 

angle-of-attack (AOA) 

deg 

4 

q 

pitch rate 

deg/sec 

5 

h 

altitude 

ft 

6 

m 

1 st flex mode 

- 

7 

m 

1 st flex mode rate 

- 

8 

72 

2 nd flex mode 

- 

9 

m 

2 nd flex mode rate 

- 

10 

73 

3 rd flex mode 

- 

11 

73 

3 rd flex mode rate 

- 


Table 1. States for Hypersonic Vehicle Model 


» 

Symbol 

Description 

Units 

1 

FER 

stoichiometrically normalized fuel equivalence ratio 

- 

2 

6e 

elevator deflection 

deg 

3 

5c 

canard deflection 

deg 


Table 2. Controls for Hypersonic Vehicle Model 


In the above model, we note that the rigid body motion impacts the flexible dynamics through the generalized 
forces. As discussed earlier, the flexible dynamics impact the rigid body motion through thrust, lift, drag, 
and moment. Nominal model parameter values for the vehicle under consideration are given in Table 3. 
Additional details about the model may be found within the following references 7-9,18-27 . 


Parameter 

Nominal Value 

Parameter 

Nominal Value 

Total Length (L) 

100 ft 

Lower forebody angle (til) 

6.2° 

Forebody Length (Li) 

47 ft 

Tail angle (r 2 ) 

14.342° 

Aftbody Length (L 2 ) 

33 ft 

Mass per unit width 

193.7107 slugs/ft 

Engine Length 

20 ft 

Weight per unit width 

6232.5 lbs/ft 

Engine inlet height h, 

4.5 ft 

Mean Elasticity Modulus 

8.6482 x 10 7 psi 

Upper forebody angle (tiu) 

3° 

Moment of Inertia I yy 

86,845 slugs ft 2 /ft 

Elevator position 

(-85, -3.5) ft 

Center of gravity 

(-55,0) ft 

Diffuser exit/inlet area ratio 

0.15 

Elevator Area 

17 ft 2 

Titanium Thickness 

9.6 in 

Nozzle exit/inlet area ratio 

6.67 

First Flex. Mode (u ni ) 

19.7 rad/s 

Second Flex. Mode (w„ 2 ) 

47.5 rad/s 

Third Flex. Mode (w„ 3 ) 

94.4 rad/s 

Flex. Mode Damping (£) 

0.02 


Table 3. Vehicle Nominal Parameter Values 


III. Elevator Surface Studies/Design 

The elevator surface is responsible for managing the unstable dynamics of the vehicle. An appropriately 
designed elevator surface will allow mission maneuvers to be executed while maintaining enough control au- 
thority to suppress un- modeled vehicle dynamics and disturbances. An overly large/poorly placed elevator 
surface will contribute unnecessary drag and add additional weight to the vehicle. The following constraints 
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must be considered when sizing/placing the elevator surface: 

Elevator Saturation. Elevator constraints such as saturation can result in performance degradation, and 
possibly instability. 45 Anti-windup and error governor methodologies may used to maintain stability at the 
cost of performance. In general it is desirable to design control surfaces such that saturation does not occur 
for standard maneuvers. Previous works have used saturation levels of [— 15°30°] and [— 20°20°] for illustra- 
tive purposes. 48,59 Instead of focusing on a particular saturation limit, this paper will attempt to show the 
tradeoffs in sizing and location for a range of deflection values. 

For the altitude pull-up maneuver, elevator usage will be largest due to one of the following scenarios: 

1. The dynamic effect of the controller executing the beginning of the maneuver. The size of the ele- 
vator deflection will depend upon (1) the magnitude of the FPA command and (2) the elevator-FPA 
bandwidth. 

2. At the lowest dynamic pressure of the maneuver. The pitching moment produced by the elevator is 
linearly proportional to the dynamic pressure. 32,49 The lOx decrease in dynamic pressure will need an 
appropriate increase in elevator deflection to maintain vehicle equilibrium. 

Rate Limitation. Rate limitations exist for the elevator due to the hydraulic actuators that manipulate 
the surface. If the BW associated with these actuators is less than an order of magnitude above the control 
BW, it should be accounted for in the control design. 37 

Separation. Elevator effectiveness at high speeds and large deflections can be over-predicted by the invis- 
cid/viscous theory due to boundary layer separation. 60 Tailing edge controls such as an elevator experience 
separation along the compression side of the surface, reducing the pressure loads. The effective force gener- 
ated by the surface are unpredictably affected by sudden changes in the pressure magnitude as separation 
regions form and dissipate. In addition, heat transfer rates can be greatly increased due to vortices within 
the separated region. 

Elevator Effectiveness. Figure 2 show the moment coefficient for the elevator vs the elevator deflection 
(AO A = 0). The figure can be utilized to determine the increase in effectiveness needed to reduced a 
particular deflection. For example to reduced a 20° deflection by 15° requires a increase in effectiveness by 
a factor of 57 % . This would need to be accomplished through some combination of an increase in elevator 
surface area and moment arm. Subsections III. A and III. A. 2 summarize the effect on the vehicle due to the 
sizing and placement of the elevator surface in greater detail. 



Figure 2. Elevator Moment Coefficient vs Deflection 


III. A. Elevator Area 

In this section we examine the impact of varying the elevator surface area. The following assumptions were 
made: 
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The elevator area is increased in steps of 1.7 ft 2 from 8.5 ft 2 to 34 ft 2 . 

The elevator location is fixed at 85 ft from the nose, 3.5 feet above the waterline. 
The weight and structural limitations of the elevator are not modeled. 

The elevator is assumed to be a flat plate; flow separation is not modeled. 
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III.A.l. Impact of Elevator Area on Static Properties (Level Flight) 

Trimmable Region. Figure 3 shows how the trimmable region depends upon elevator area. 


Envelope Variations with Elevator Area 



From Figure 3, the following changes in the trimmable region are apparent: 

• The trimmable region shrinks with decreasing elevator area 

• The pinch point decreases in Mach and altitude with decreasing elevator area 
Trim AO A. Figure 4 shows how the AO A depends upon elevator area. We see that 

• Trim AOA is almost independent of the elevator area. 


AOA vs. Elevator Area, h=100 kft 


AOA vs. Elevator Area, M=8 
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Figure 4. Trim AOA vs Elevator Area 


Trim Elevator. Figure 5 shows how trim elevator deflection depends upon elevator area. We observe that 

• Trim elevator deflection decreases with increasing elevator area. 

• As the elevator area increases, the deflection is weakly dependent on the operating point. 

Lift-to-Drag ratio. Figure 6 shows how trim L/D ratio depends upon elevator area. We observe the following: 
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Elevator Inviscid Drag (Ibf) • Lift to Drag ratio Elevator deflection (deg) 


Elevator vs. Elevator Area, h=1 00 kft 


Elevator vs. Elevator Area, M=8 



Figure 5. Trim Elevator vs Elevator Area 



At lower altitudes, the L/D ratio is concave down with respect to (w.r.t) the elevator area. Figure 5 
shows that the elevator deflection decreases with increasing elevator area. The inviscid drag decreases 
with decreasing elevator deflection (see figure 7). However, the elevator’s viscous drag increases. Hence 
the total elevator drag follows a non-monotonic (concave up) trend w.r.t. elevator area (figure 8). 



Trim FER. Figure 9 shows how trim FER depends upon elevator area. From Figure 9, one observes that 
the: 

• Trim FER is non-monotonic (concave up) with increasing elevator area at low altitudes (since total 
drag is concave up w.r.t. elevator area at low altitudes). 
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Total Elevator Drag vs. Elevator Area, M=8 



Figure 8. Elevator Drag vs Elevator Area 


FER vs. Elevator Area, h=100 kft 




10 15 20 25 30 35 

Elevator Area (ft 2 ) 


vs Elevator Area 


III. A. 2. Impact of Elevator Area on Dynamic Properties 

Instability. Figure 10 shows how the instability varies with the elevator area. From Figure 10, one observes 


RHP Pole vs. Elevator Area, h=100 kft RHP Pole vs. Elevator Area, M= 8 



Figure 10. Right Half Plane Pole vs Elevator Area 


that the: 

• The instability decreases with increasing elevator area - a larger elevator moves the aerodynamic center 
(ac) rearward, and closer to the center of gravity (eg). As a result, the instability decreases. 

• As the elevator area increases, the dependence of the RHP pole on the operating Mach decreases. 

RHP Zero. Figure 11 shows how the RHP zero depends varies with the elevator area. From Figure 11, one 
observes that the: 

• RHP Zero remains relatively constant with increasing elevator area 


12 of 23 


American Institute of Aeronautics and Astronautics 



Elevator Area (ft 2 ) Elevator Area (ft 2 ) 

Figure 11. Right Half Plane Zero vs Elevator Area 


III.B. Elevator Location 

In this section we examine the impact of moving the elevator rearward w.r.t the the vehicle nose. The 
following assumptions were made: 

• The horizontal elevator location is varied from 70 to 90 ft from the vehicle nose. 

• The vertical elevator location is fixed at 2.5 ft above the waterline. 

• The elevator area is fixed at 17 ft 2 

• The weight and structural limitations of the elevator are not modeled. 

• Elevator is assumed to be a flat plate; flow separation is not modeled. 

III.B. 1. Impact of Elevator Location on Static Properties 

Trimmable Region. Figure 12 shows how the trimmable region depends upon elevator location. Specifically, 
we see that the 

• Trimmable region expands with as the elevator is moved rearward. 

• Pinch point increases with Mach and altitude as the elevator is moved rearward. 

Trim AOA. Figure 13 shows how the AOA depends upon elevator location. We see that 

• Trim AOA increases as the elevator is moved rearward. 

Trim Elevator. Figure 14 shows the variation in trim elevator deflection with a rearward elevator. We 
observe that the 

• Trim elevator deflection decreases as the elevator is moved reward - the elevator moment arm increases, 
and hence less elevator deflection to provide the same moment. 

• As the elevator is moved rearward, the dependence of the trim elevator deflection on Mach decreases. 

Lift-to-Drag ratio. Figure 15 shows how trim L/D ratio depends upon elevator location. We observe the 
following: 

• Trim L/D ratio increases as the elevator is moved reward - less elevator deflection for a fixed elevator 
area results in less drag and higher trim L/D. 

FER. From figure 16, we find that the 

• Trim FER decreases as the elevator is moved reward - less elevator deflection for a fixed elevator area 
results in less drag and less trim FER. 
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Mach 

Figure 12. Trimmable Region Dependence on Elevator Location 


AOA vs. Elevator Location, h=1 00 kft 


AOA vs. Elevator Location, M=8 



Figure 13. Trim AOA vs Elevator Location 


Elevator vs. Elevator Location, h=100 kft 


Elevator vs. Elevator Location, M=8 



Figure 14. Trim Elevator vs Elevator Location 


III.B.2. Impact of Elevator Location on Dynamic Properties 

Instability. Figure 17 shows how the vehicle instability depends upon the elevator location. It is seen that 
the: 
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RHP Pole • FER Lift to Drag ratio 


L/D vs. Elevator Location, h=100 kft L/D vs. Elevator Location, M=8 



Figure 15. Trim L/D vs Elevator Area 


FER vs. Elevator Location, h-100 kft FER vs. Elevator Location, M=8 



instability decreases linearly as the elevator is moved rearward - the ac moves rearward with the elevator 
(and closer to the eg). As a result, stability imroves. 


RHP Pole vs. Elevator Location, h=1 00 kft RHP Pole vs. Elevator Location, M= 8 
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RHP Zero. Figure 18 shows how the vehicle instability depends upon the elevator location. From 



Figure 18, one observes that the: 

• RHP zero increases linearly as the elevator is moved rearward, due to the increasing moment arm of 
the elevator. 

III.C. Summary - Aerodynamic control surface studies 

Elevator Area Increase. Table 4 summarizes trends for increasing elevator areas. 61 As the elevator area is 
increased, we observe specific monotonic trends and tradeoffs that result in the following PROS and CONS: 

• PROS: trim AOA nearly constant, trim lift remains nearly constant, trim elevator decreases, RHP pole 
decreases, RHP zero nearly constant. 

We also observe the following more complex (non-monotonic) behavior: 

• trim lift-to-drag increases (max depends on Mach, altitude) and then decreases; 

• trim drag, FER and fuel usage decreases (min depends on Mach, altitude) and then increases. 


Table 4. Trends for Increasing Elevator Area 


Property 

Pro 

Con 

Trim Lift 

Almost Constant 

Almost Constant 

Trim Drag (minimizer) 

t with Alt. (fixed Mach) 

4_ with Mach (fixed Alt.) 

Trim L/D (maximizer) 

t with Alt. (fixed Mach) 

with Mach (fixed Alt.) 

Trim AOA 


Almost Constant 

Trim Elevator Decreases 



Trim FER (minimizer) 

t with Alt. (fixed Mach) 

4. with Mach (fixed Alt.) 

Trim Fuel Rate (maximizer) 

t with Alt. (fixed Mach) 

4. with Mach (fixed Alt.) 

RHP Pole 

Decreases 


RHP Zero 


Almost Constant 

RHP Z/P Ratio 

Increases 



Rearward, Elevator Shift. Table 5 summarizes trends for varying the elevator location. 61 As the elevator is 
moved rearward, we observe specific monotonic trends and tradeoffs that result in the following PROS and 
CONS: 

• PROS: trim L/D, RHP zero and RHP zero-pole ratio increase monotonically. Trim elevator, FER, 
RHP pole decrease monotonically. 
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• CONS: trim AOA, fuel rate increase monotonically. 
Non-monotonic patterns are not observed in this case. 


Table 5. Trends for Rearward Elevator Shift 


Property 

Pro 

Con 

Trim Lift 

Almost Constant 

Almost Constant 

Trim Drag 

Decreases monotonically 


Trim L/D 
Trim AOA 

Increases monotonically 

Increases monotonically 

Trim Elevator 

Decreases monotonically 


Trim FER 

Decrease monotonically 


Trim Fuel Rate 
RHP Pole 

Decreases monotonically 

Increases monotonically 

RHP Zero 

Increases monotonically 


RHP Z/P Ratio 

Increases monotonically 



IV. Control System 

A classic decentralized inner-outer loop control system architecture was used to illustrate control design 
issues. Such an architecture was examined within [43]. It can be visualized as shown in Figure 19. 



Figure 19. Inner Outer Loop Control System 


Outputs: y p = [ v 7 ] T Controls: u = [ FER 5 e ] T State Feedback: x r = [9 0] T 


Overview of Nominal Control System Architecture. Within Figure 19, we have an inner loop con- 


troller: 



0 


0 


- 9A 

s 2 -|-bi s-\-b 0 

-9i^ 

s 2 -\-bis-\-b 0 


s 2 +ais+a 0 

s 2 +ais+a 0 

_ 


( 9 ) 


and an outer loop controller: 


K 0 (s) 


9b(s+Zb ) 
s 


10 

s+10 


1 3 


0 


g-f(s+z y 


20 


s+20 


(10) 


The inner-outer loop structure is now described. 
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1. Velocity Control Via Single Loop PI. A PI controller (with roll off) is used to control velocity. Two 
parameters are associated with velocity control: ( gb , Zb ). A simple anti- windup method (not discussed 
in the paper) is used to address the state-dependent thermal choking nonlinearity discussed earlier. 

2. FPA Control Via Inner-Outer Loop. A PD (proportional plus derivative with notch) inner-loop con- 
troller is used on pitch to stabilize the vehicle’s pitch dynamics and make the modified dynamics look 
friendly to the outer-loop FPA PI (proportional plus integral) controller. The inner loop’s notch pro- 
vides lead to assist with the first flexible mode’s complex zero-pole pair (near 21 rad/sec). As might 
be expected with this inner-outer loop structure, the inner loop is generally much faster than the outer 
loop. Its bandwidth must be large enough to stabilize the vehicle and follow outer-loop commands, but 
not too high so that the flexible dynamics are overly excited. (It should be noted that any practical 
“real-world” implementation of the PD controller should include additional roll off in order to attenuate 
high frequency noise rq resulting from 9 and theta measurements.) 

A PI (proportional plus integral with roll off) outer-loop controller is used for FPA. The outer loop is 
generally much slower than the inner loop. The maximum achievable outer loop bandwidth is limited 
by the vehicle’s flexible dynamics as well as the right half plane zero associated with the elevator-FPA 
map. 

Four parameters are associated with FPA control: (g t , Zi ) for the inner-loop PD controller, g 7 , z 7 ) for 
the outer-loop PI controller. 

Command pre-filters are also included: IT), = -f*- - for velocity reference commands, Wfpa = 2Z for FPA 

s~r Zfj s~r z~y 

reference commands. (These reference command pre-filters are not shown within Figure 19.) 

Additional information about controller parameter selection can be found within [43]. Nominal control 
system parameter values are given in Table 6. 


9i 

Zi 

9i 


9b 

Zb 

o 0 

a\ 

bo 

in 

1 

4.5 

21.23 

0.404 

2.2437 

0.0278 

30 

900 

15 

250 


Table 6. Control System Parameters for Nominal Hypersonic Vehicle (Mach 8, 85kft, Level Flight) 

Overview of Nominal Guidance System Architecture. 

Within this paper, we consider a constant q = 2076 psf flight profile followed by a pull-up maneuver. A 
guidance system is necessary to issue appropriate FPA commands to maintain a constant dynamic pressure. 
The complete guidance and flight control systems may be visualized as shown below in Figure 20. The 
guidance system consists of a PI feedback structure (with roll off): 

K guidance (s) = \_Zl1 ( 11 ) 

S J L s -I- z 1 _ 

where k = 1/7846, z = 15.7. z 1 is determined from the outer-loop FPA controller - it is roughly selected to 
be the FPA outer-loop’s PI controller’s zero or larger. The second term in brackets is the FPA command 
pre- filter. 

• A simple PI guidance system (which processes dynamic pressure error = (72076 psf — factual) is used 
to generates FPA guidance commands until Mach 8, 85 kft; 

• A step velocity reference command is used to take the vehicle along q = 2076 psf from Mach 5.7 at 70 
kft to Mach 8 at 85 kft. 

• Once Mach 8, 85 kft is reached, the pull-up maneuver is executed; FPA commands of increasing 
magnitude will be compared. 

For a more detailed discussion on gain scheduling of the nominal design, see [43]. Within [62] the author 
also examines gain scheduling issues for the nominal vehicle model. 
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rq = 2076 psf 





^ guidance \& ) 


e q 



Figure 20. Inner-Outer Loop Architecture and Guidance System 


IV. A. Nonlinear Simulation 

Within this section we consider an altitude pull-up maneuver initiated at Mach 8, 85 kft. FPA commands 
magnitudes of 3°, 5°, and 7° have been issued to execute the maneuver. A nominal elevator area of 17 ft 2 
has been chosen. The elevator has been place at 85 ft back from the nose of the vehicle and 3 ft above the 
water line. The CG of the vehicle is located at approximately 55 ft from the nose of the vehicle along the 
water line. Important results have been summarized within Table 7. 


OUTPUTS 




Figure 21. Pull-up Trajectories: Mach and FPA 
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Figure 21 shows the Mach number and commanded FPA commands for the 3 trajectories. Mach is com- 
manded to remain at 8, but cannot be maintained there due to thermal choking limits within the engine 28,48 
as seen in Figure 22. 


CONTROLS 


FER tc 



















O' 

0 20 40 60 80 100 120 

t (sec) 



Figure 22. Pull-up Trajectories:Controls 


For 3° and 5° FPA commands, a peak elevator usage of 20° occurs when the vehicle reaches q = 200 psf. If 
the saturation limit for the elevator was 20°, and a 5° safety margin was desired for disturbance rejection, 
then the elevator effectiveness would need to be increased by roughly 57 % as detailed in Section III. For 
the 7° (and larger) command, peak elevator usage occurs at the beginning of the pull up maneuver. Peak 
deflection at this location can be reduced through a combination of improved elevator effectiveness and 
decreased closed loop bandwidth. 


STATES 



140 



0 20 40 60 80 100 120 

t (sec) 


Figure 23. Pull-up Trajectories: AO A and Altitude 
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Table 7. Trajectory Properties 



V. Summary and Directions for Future Work 

Summary and Conclusion. Elevator surface studies have been conducted w.r.t. size and placement. The 
impact on the vehicles trimmable envelope (region in altitude-Macli space), static (equilibrium) conditions 
at trim, and dynamic properties as they relate to control system design was addressed for level flight. These 
studies demonstrates the tradeoffs that can be expected if elevator deflection needs to be decreased due to 
saturation constraints. Elevator usage has been shown for a family of pull-up trajectories. 

Directions for Future Work. Future work will examine the impact of reduced elevator effectiveness due 
to phenomena such as separation. Anti-windup and error governor will be implemented within the control 
law to examine their effectiveness during pull-up maneuvers. 

Nomenclature 

true air speed 
angle of attack (AO A) 

flight path angle (FPA); not to be confused with specific heat ratio for air 
pitch angle 
pitch rate 
altitude 

generalized elastic coordinates (i = 1,2,3) 
generalized elastic rates (« = 1, 2, 3) 
elevator deflection 

normalized fuel equivalence ratio (£) 
flow deflection angle 
Lower Forebody Wedge Angle (6.2 deg) 
density of air (standard atmosphere) 
dynamic pressure (^pv 2 ) 
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